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Two electronic channels competing to screen a single impurity spin, as in the two-channel Kondo model, are 
expected to generate a ground state with nontrivial entanglement structure. We exploit a spin-chain representa¬ 
tion of the two-channel Kondo model to probe the ground-state block entropy, negativity, tangle, and Schmidt 
gap, using a density matrix renormalization group approach. In the presence of symmetric coupling to the 
two channels we confirm field-theory predictions for the boundary entropy difference, \n(gyy jgjf) = ln(2)/2, 
between the ultraviolet and infrared limits and the leading ln(x:)/x impurity correction to the block entropy. 

The impurity entanglement, 5imp, is shown to scale with the characteristic length ficK- We show that both the 
Schmidt gap and the entanglement of the impurity with one of the channels - as measured by the negativity- 
faithfully serve as order parameters for the impurity quantum phase transition appearing as a function of channel 
asymmetry, allowing for explicit determination of critical exponents, vx2 and fx 0.2. Remarkably, we find the 
emergence of tripartite entanglement only in the vicinity of the critical channel-symmetric point. 


Introduction.- The Kondo effect is one of the most intrigu¬ 
ing effects in quantum many-body physics. At low tempera¬ 
tures, a localized magnetic impurity is screened by the con¬ 
duction electrons leading to the formation of many-body en¬ 
tanglement. A generalization of the Kondo model was intro¬ 
duced by Nozieres and Blandin m, where another channel 
of electrons is also coupled to the impurity. This is the well- 
known two-channel Kondo (2CK) model, for which various 
results were obtained using Bethe ansatz conformal 

field theory ||5] ID (CFT), bosonization ||7}j9l and entangle¬ 
ment of formation Qo). This model is very different from the 
one-channel Kondo (ICK) model as the two channels compete 
to screen the spin-1 /2 impurity, leading to an “overscreened” 
residual spin interacting with the electrons 0. This leads 
to non-trivial properties including a residual zero-temperature 
impurity entropy and a logarithmic behavior of magnetic sus¬ 
ceptibility and specific heat. However, channel symmetry is 
crucial; even the smallest asymmetry leads to screening of the 
impurity by the channel with the stronger coupling 0, and as 
the channel asymmetry is varied, an impurity quantum phase 
transition (IQPT) occurs at the symmetric point, correspond¬ 
ing to the 2CK model. 

Intensive research has been carried out to investigate the 
thermodynamics and the transport properties of the 2CK 
model ifTHSl l5H9l rTTH25]l . Experimentally, signatures of the 
2CK model have been observed in mesoscopic structures M- 
|29l . Still, the real-space entanglement structure and the im¬ 
prints of the two distinct length scales ^ 2 Ck ~ m/7’2ck and 
~ m/T* with u the spin velocity- implied by the known 
crossover energy scales 7’2 ck (2CK temperature) and T* (crit¬ 
ical crossover in the channel-asymmetric case) 0[23 - have 
not yet been unraveled. A way forward is to use a spin-chain 


representation of the 2CK model imiia , which allows for ef¬ 
ficient Density Matrix Renormalization Group (DMRG) com¬ 
putations USSEl {m - 100 - 1024 states kept) to uncover the 
ground state entanglement properties. 

In this letter, we show how the implementation of this 
scheme allows for a detailed study of the entanglement in 
the 2CK ground state and the IQPT between the two channel- 
asymmetric ICK phases. Specifically, we present results for 
the impurity entanglement entropy ll^ [3TI| . the negativity 
IITSIIMII . the Schmidt gap llJTlI^ . and the tripartite entangle¬ 
ment ll^l40ll . At the channel-symmetric 2CK point we show 
that ^ 2 Ck can be interpreted as a dynamically generated cut¬ 
off length by demonstrating scaling of the impurity entangle¬ 
ment entropy. A detailed analysis allows us to extract the two- 
channel boundary entropy difference \t\{giiv!gm) - ln(2)/2, 
between the ultraviolet and infrared limits ED, as well as the 
leading correction ln(x)/x, for block sizes x » ^ 2 Ck ED- In 
addition, we show that the negativity and the Schmidt gap act 
as order parameters for the IQPT, enabling us to predict, via 
finite-size scaling, the pertinent critical exponents. Finally, we 
compute the tangle Il3^l40ll and show that tripartite entangle¬ 
ment emerges only in the vicinity of the critical point. 

Spin-chain representation.- We consider two open Heisen¬ 
berg chains coupled to a single spin-1/2 impurity as shown 
in Fig. [2 a). The open chain Hamiltonian is given by 
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Neumann entropy. By doing so, we provide a direct ’’quan¬ 
tum probe” of the boundary entropy predicted by CFT ED, 
with no reference to the thermodynamic entropy. 

Let’s first consider the N ^ co limit. When J' - 1, we 
simply have a uniform periodic chain with region A consisting 
of 2r H- 1 sites. Then, using the fact that the central charge 
c = 1, the entanglement entropy for region A of a periodic 
chain is predicted to be, from CFT Il46l 

5a(/= l,x,A^) = ^ln(2rH-1)-Hsi (3) 


FIG. 1. (Color online) (a) Kondo spin chain with a spin-1/2 impurity cou¬ 
pled to its left and right channels by TJ' and J', respectively. For F = 1 the 
impurity is screened by both channels representing the 2CK model while for 
r 1 ICK physics emerges, (b) The impurity entropy Simp is computed as 
the difference between the entropy of region A with and without the impurity, 
(c) Partitioning of the system for computing the Schmidt gap. 

the impurity spin and the spin at site I in channel m, respec¬ 
tively, and Nm is the number of spins in chain m making the 
total number of spins N - Nl + Nr + \. We choose the 
nearest-neighbor coupling Ji to be unity and the next-nearest- 
neighbor coupling J 2 - J 2 (with - 02A\2J\) so as to 
remove marginal coupling effects Il34ll43l . In this work, we 
set the Kondo coupling 7^ = F7' and 7^ = 7', with F = 7^/7^ 
keeping 7,'„ < 1. The Hamiltonian has been introduced in 
Ref. ifTTl fT2ll as a representation of the spin sector of the 2CK 
model when F = 1. For further justifications, see the Supple¬ 
mental Material. For any F + 1 ICK physics emerges. For the 
case of F = 1, we also use a periodic chain, as shown on the 
left side of the Fig.[TJb), by adding the following terms: 

HpBC - Hqbc + ■ O’** + J'JiO'b ■ O"^ 

+ 7 , (<r^. (2) 

Again N = Nr + Nr + I, and at 7' = 1 we obtain a uniform pe¬ 
riodic chain which presents significant advantages Il44ll . In the 
limit of A ^ 00 the two boundary conditions are equivalent. 
For Hqbc the parity of Nl - Nr is crucial m but here we 
only study Nl - Nr odd, however, for Hrbc it is the parity of 
N that matters m and we only study N even (Nl - Nr ± 1) 
which makes the parity effects compatible for Hqbc and Hrbc- 
Impurity entanglement entropy.- We first study the channel- 
symmetric case, F = 1, with ^ 2 Ck being the only relevant 
length scale in the problem. We consider the von Neumann 
entropy, Sa(J',x, N) = -Trp^t log pA with pA the reduced den¬ 
sity matrix of a region A which includes the impurity spin and 
X spins on either side of it. N is the total number of spins in the 
system, including the impurity. We consider an even periodic 
system, using Hrbc as shown in Fig.[TJb). This boundary con¬ 
dition should not affect our results as long as r A/2 El- 
Similar to the single-channel case ll^ l3T1l the entanglement 
entropy behaves very differently in the two limits x ^ 2 Ck 
and X » ^ 2 CK, with ^ 2 Ck ~ growing exponentially as 
7' ^ 0 (for some constant d). In what follows we shall show 
how to pinpoint the impurity contribution. Simp, to the von 


for a non-universal constant si. For finite but large A even, we 
expect the limit of 7' —> 0^, x N (which is different from 
the case where the impurity is absent) to give: 

5A(7'^0+,r,A) = 5A(r,A-l)H-ln2 (4) 

where S a{x,N-1) represents the entropy of region A when the 
impurity is absent but the region still consists of x spins from 
each channel (so the total length is A-1) as shown on the right 
side of Fig. [^b). The additional In 2 entanglement entropy in 
Eq. Q is the impurity contribution and can be understood by 
observing that a spin chain with an even number of sites has a 
spin zero ground state for any 7' > 0 no matter how small. In 
a valence bond picture of the A even ground state there will 
always be an (impurity) valence bond (IVB) connecting the 
impurity spin to another spin in the system, although the IVB 
becomes very long in the small 7' limit ll^ l3TI . Intuitively, 
this long IVB adds an extra In 2 to SaW 0^,x, A). The 
interesting case of A odd will be considered elsewhere l45l . 

In the absence of an impurity, as long as x « A/2, the en¬ 
tropy of region A is the sum of the entropy of two equal blocks 
at either end of an open chain as shown in the right part of 
Fig.[TJb). In this case the open boundaries induce an alternat¬ 
ing term in the entanglement entropy BtII and we therefore 
only focus on the uniform part, S “ finding Ebl l48l 

5“(x,A-l) = 2[iln(2x) + y+lng], (5) 

where sj is the same non-universal constant appearing in 
Eq. ([^ and In g is a universal term arising from a non-integer 
“ground-state degeneracy”, g ED- 

The difference between the two entropies of the two ex¬ 
treme regimes will be 

Sa(J' - l,x,A)-5^(7' —> 0^,x,A) = -21ng-ln2-i-(9(l/x). 

( 6 ) 

Using the mapping of the spin-chain system onto the 2CK 
model, we associate 7' —> 0^ with the weak coupling ultra¬ 
violet fixed point and 7' —» 1 with the infrared fixed point. 
Hence we expect 

5^(7' = l,x,A)-5“(7'^0+,x,A) = lng 2 R-lngj/v, (7) 

where In guv and Ing/s are the boundary entropies for the ul¬ 
traviolet and infrared fixed points. Hence, it follows that the 
degeneracies of the 2CK model and the open chain must be 
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related as guvIgiR - 2g^- While guv - 2, corresponding 
to the decoupled impurity spin, g/g has the non-trivial value 
of V 2 . On the other hand, g was predicted, using field the¬ 
ory arguments for the open spin chain to have the value 2 
|[TT][T2]|49] validating the relation = 2g^ = V 2 . This 

constitutes a highly non-trivial check of the spin-chain repre¬ 
sentation of the 2CK model. We confirm the result g - 2 
by extracting ii from DMRG results for the entanglement en¬ 
tropy for an even periodic chain, finding si = 0.743743. We 
then fit 5 ^ for a single open chain of length N to the finite N 
generalisation of Eq. 0 


5“(x,iV) = 2 


- ln[(2A^/7r) sin(;7rx/A^)] - 1 —-H Ing 
6 2 


a, 


+ — [2 H- 7 r(l - 2xlN) cot{nx/N)]. 


( 8 ) 


Here the last term is a correction, behaving as ajx in the N —> 
00 limit, calculated in OOl [311 l45ll where a is a non-universal 
parameter. is extracted using a 7-point formula li^ [311 
[ 45 I . With Si known, this then determines Ing = -0.17328, in 

excellent agreement with ln(2“i) = -0.1732867_ 

We now show that \vi{giiv!gm) enters as part of the limit¬ 
ing behavior of the impurity entanglement entropy allowing 
us to numerically estimate this boundary entropy difference. 
We begin by considering the behaviour of 5 .4 for intermediate 
values of J'. Most notably, an alternating term appear in 5 a 
for any 7' 4 1 1[45) . Hence, by subtracting off the entropy with 
the impurity absent IISOIIMI . as shown in Fig.[ 2 b), we define 
the impurity entanglement entropy using the uniform part as 


5i„,p(7',x,fV) = 5"(/,x,fV) - 5“(x,fV- 1). (9) 


The hallmark feature of the characteristic length ^ 2 Ck ~ 
M/T 2 CK is that Simp is a universal scaling function of the two 
variables x/N and x/f 2 CK- Again, the parity of N also plays 
a crucial role 1451 but here we only focus on N even. If we 
fix xjN - 1/10, Simp should then be a function of the single 
variable xjfxcK- However, as evident from Eq. ([^ the term 
proportional to a in S‘^(x,N - 1) gives rise to corrections to 
scaling disappearing as 77 ^ 00 with x/N and x/f 2 CK held 
fixed. Eor clarity, we subtract these corrections from S imp ob¬ 
taining Sjl^. In Eig. I^a) we demonstrate the scaling by col¬ 
lapsing data for many values of J' and N with fixed x/N onto a 
single curve by appropriately selecting ^ 2 Ck( 7 ')- The expected 
^ 2 CK ~ behavior is also confirmed (inset of Eig.|^a)). We 
see that an excellent data collapse appears for a range of J' and 
the data approaches fairly closely to ln(2)/2 = 0.3465 at large 
x/ficK ■ This limit corresponds to 7' ^ 1 and using Eq. 0 
and 0 we have S imp(7' ^ 1) = ln(2)-ln(g£/y/g/«) = ln(2)/2 
so we can conclude l[44l : 

^^iguv/gm) - ln( 2 )/ 2 , guv/grn - V 2 ( 10 ) 


providing a firm confirmation of the CET predictions. 

Eor X » ^ 2 Ck at 77 ^ cm we are close to the infrared fixed 
point. The leading irrelevant operator has dimension 3/2 a 
and is expected to lead to corrections to the leading ln( 2)/2 



FIG. 2. (Color online) (a) Scaling of for fixed x/N =1/10 

(A^ even). At J' = 0.9, ^2Ck(^ 0 is arbitrarily fixed at 0.07747 to coincide with 
the estimate from panel b. Inset: ^2 Ck( 2') as a function of 1/7'. (b) DMRG 
results for 7' = 0.9,A = 800). For^2C^(7') « x « A/2 can 

be fit to the form A ln(jr/^ 2 CA:)/(^/^ 2 CA') + B (red line) with ^ 2 Ck{J' = 0-9) = 
0.07747, A = 0.69 and B = 0.34 ~ ln(2)/2 significantly better than to ~ l/x 
(green line). Inset: Convergence to the limiting form at x A/2 with A. 


behavior of 5 imp that in second-order perturbation theory are 
of the form imp oc ln(x)/x El, valid in the regime ^ 2 Ck ^ 
X « 77/2. Numerically we can confirm this by studying S 
for 7' ~ 1 where ^ 2 Ck is small. This is shown in Fig. |^b) 
for 7' = 0.9 where a fit to ln(x)/x correction is statistically 
superior to a simpler l/x form over a significant range of x. 

Negativity as an order parameter.- Several entanglement 
measures have been used to detect quantum phase transitions 
ElESlISoHSl. Here, we propose the negativity |[T5l |36l as 
an order parameter for the IQPT, with E as control parameter. 
For any bipartite density matrix pab the negativity, as an en¬ 
tanglement measure, is defined as Ea.b = “1 + Itlkl, where 
Tjk’s are the eigenvalues of the matrix where Ta stands for 
partial transposition with respect to subsystem A EU. In this 
section, and through the remainder of the paper we use Hqbc, 
with Ni - Nr odd. In Fig.j^a) we plot the negativity between 
impurity and right channel, Eq r, versus E. It is expected that 
the ground state is overscreened only at E = 1 where the impu¬ 
rity is entangled with both channels. For any E 4^ 1 in the ther¬ 
modynamic limit, the impurity is screened only by the channel 
with the strongest coupling to the impurity, resulting in a fully 
screened ICK phase. Indeed, the behavior of the negativity 
is consistent; Eo.r goes from 1 to 0 around the critical point. 
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FIG. 3. (Color online) (a) Negativity between the impurity and the right 
channel (i.e. Eo r) versus T for N = 403 and J' = 0.4. (b) Derivative of Eo r 
with respect to T for different system sizes, (c) Finite size scaling of Eor. (d) 
Finite-size scaling of the Schmidt gap. 


The thermodynamic behavior can be explored by studying the 
derivative of the negativity with respect to T, namely 
shown in Fig. [^b). As the figure shows, the derivative dips 
at the critical point with the dip sharpening as N increases. 
This suggests that, as A —> oo, diverges at the critical 
point, implying that the 2CK ground state is destroyed and 
ICK physics is emerging. 

The interpretation of the negativity as an order parameter 
can be justified by a finite-size scaling analysis Il55l . An order 
parameter scales as IF - 1in the vicinity of the critical point 
and the correlation length as |F-11 where fi and v are critical 
exponents. The role of a correlation length is here taken by the 
critical crossover scale at which the renormalization-group 
flow of the channel-asymmetric model crosses over from the 
unstable overscreened fixed point to the fully screened Kondo 
fixed point siia. Finite-size scaling ll55l implies that 

1|A‘^’'), (11) 

with F a scaling function. In Fig. [^c), we plot A^^’'£’o,r as a 
function of(F- DA'-'’'. When v = 2 + 0.05 and= 0.2+0.02, 
curves for different A collapse to a single curve. The value of 
v^2 matches CFT ||3 and bosonization results |T9l , verifying 
that the negativity behaves as an order parameter. Here v = 
Ijd, with d the scaling dimension of the relevant operator that 
appears in the Hamiltonian when parity symmetry is broken. 

Schmidt gap. - Another key quantity, related to the entangle¬ 
ment spectrum, is the Schmidt gap Ay. Given a bipartitioning 
of the system, it is defined by Ay = Ti - A 2 , where Ai > A 2 
are the two largest eigenvalues of the reduced density matrix 
of any of the two subsystems. It was recently shown that the 
Schmidt gap can serve as an order parameter across quantum 
phase transitions llJTlI^ . For the 2CK model close to F = 1, 
and choosing a bipartition as shown in Fig. [^c) for two com- 




FIG. 4. (Color online) (a) The tripartite entanglement indicator r versus F. 
(b) Entanglement between the two channels versus F. In both panels /' = 0.4. 


plementary left and right blocks, the Schmidt gap is found to 
obey finite-size scaling with the same critical exponents as the 
negativity. Fig.|^d) shows the Schmidt gap data collapse for 
three different system sizes, confirming it as an alternative or¬ 
der parameter to the negativity in the 2CK model. 

Tripartite entanglement.- Changing from ICK to 2CK 
physics changes the entanglement structure fundamentally. 
Inspired by tangle 091 and its generalization for negativ¬ 
ity iQl, as tripartite entanglement measures for qubits, we 
introduce a tripartite entanglement indicator as 

T = (tto + + ^r)/3 (12) 


in which 


Z7’2 Z7’2 Z7’2 Z7’2 Z7’2 Z7’2 

^0 - “ ^0,L ~ ^O.R’ “ ^m,0m ~ ^m,0 ~ 


where m - L,R and m — R,L represent opposite channels, 
Eq.lr = 1 is the negativity of the impurity with the rest of 
the system, Fo.m = £^m,o is the negativity between the impu¬ 
rity and channel m, and E^.m is the negativity between the 
two channels. For systems with odd length leads each chan¬ 
nel effectively behaves like a spin-1/2 system and our tripartite 
entanglement indicator r becomes a natural generalization of 
the tangle defined for three qubits OOl . In Fig. Qa) we plot 
T versus F for systems with odd length leads, t clearly peaks 
at the critical point with the peak becoming more pronounced 
with increasing length, suggesting its divergence with A. The 
emergence of tripartite entanglement is therefore related to the 
overscreening at the critical point where the two channels be¬ 
come highly entangled. In Fig.j^b), we plot the negativity be¬ 
tween the two channels, Eir, versus F. As the figure shows, 
Ei^k is maximal at F = 1, likely diverging with A. 

Conclusions.- Employing high-precision DMRG computa¬ 
tions, we have studied the ground state entanglement of the 
2CK model, allowing us to uncover the fractional ground state 
degeneracy predicted by CFT. The existence of the character¬ 
istic length scale ^ 2 ck is established through a scaling anal¬ 
ysis of Simp. The IQPT appearing as a function of channel- 
asymmetry and its exponents is detected using both the nega¬ 
tivity and the Schmidt gap as order parameters. Furthermore, 
the tangle is used to show that tripartite entanglement emerges 
only in the vicinity of the critical point. 
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Supplemental Material 

1. Review of negativity as an entanglement measure 

In a bipartite system AB, quantification of entanglement between the two subsystems has been intensively studied in the last 
decade. When the density matrix of the overall system p^g is pure then the von Neumann entropy of the subsystem, namely 
S ip a) - -Tr{pA logP a) where pA is the reduced density matrix of the subsystem A, is a unique measure of entanglement and 
all other measures are monotonic functions of the von Neumann entropy. At variance, when the overall state pab is mixed the 
von Neumann entropy fails to quantify the entanglement between the two subsystems. In such cases one can use negativity as a 
pertinent measure of entanglement, dehned as lIlSl 1^ 

E^-l+Yj\^k\^2Y,\Vk\, (SI) 

k t]t<0 

where the rj^s are the eigenvalues of the matrix pj"^ (or pjp in which Ta (or Tg) stands for partial transpose with respect to 
the subsystem A (or B). If the density matrix pab is separable, then both pj"^ and p^^ remain positive and thus tju > 0 for all 
A:’s which results in zero negativity in Eq. In contrast, if the overall state pab is entangled then some of the t/^’s become 
negative and thus the negativity in Eq. ( |Sl| l becomes nonzero. 

Negativity as a measure of entanglement, applicable for both pure and mixed states, is an entanglement monotone which 
means that it does not increase under local operations. Eurthermore, negativity is a legitimate quantum mechanical observable 
in the sense that it is associated with a Hermitian operator as 

O - 2 ^ {\m){rik\f^ , (S2) 

tik<0 


where |77jt)’s are the eigenvectors of the matrix p^^. So, one can easily show that E - Tr(pAB(9). 


2. Spin-chain representation of the 2CK model 

The spin-chain representation of the spin sector of the two-channel Kondo model was first introduced in Ref. nulla. To 
further justify this representation, we here provide a conformal held theory analysis. 

By folding the system in half, creating an open boundary at the impurity site I - 0, the left and right parts of the spin chain 
come to dehne two ’’channels” on the interval [0, Aa/2], where a is the lattice spacing. Introducing slowly varying left/right- 
moving SU(2)i spin currents together with two SU(2)i Wess-Zumino-Witten matrix helds g"' Ea, the representation 

^ + (-l)'Tr(g('”>(fl0o-), (S3) 

can then be used in the continuum limit to map the folded T = 1 system onto the spin sector of the two-channel Kondo model 
iniaa- Since the two channels (a.k.a. spin chains) both couple to the impurity, the expected chiral SU(2)i(8> SU(2)i symmetry 
of the critical theory in absence of the impurity gets replaced by SU( 2)2 ® Z 2 , where the Ising symmetry group Z 2 encodes the 
presence of the two channels. This conformal embedding is different from the one used for the original two-channel Kondo 
model 0, and reflects the fact that the two underlying bulk theories are different; In the two-channel Kondo model one is 
dealing with non-interacting electrons with two orbital channels, whereas in its spin-chain representation the impurity host 
is that of a one-dimensional half-hlled Hubbard model (with charge gapped out) and with only a single orbital. The Ising 
Z 2 sector comes into play only off the two-channel Kondo critical point, and then contributes a leading scaling operator of 
the same dimension as that from the ’’flavor” sector in the conformal embedding used for the two-channel Kondo model 0. 
This testihes to the consistency of the spin chain representation. Note that the total central charge implied by the spin-chain 
representation on a half-line, c - csu( 2 ) + cising = 3/2-i-l/2 = 2, becomes halved, c - 1, when doubling the spatial degrees of 
freedom by unfolding the system back to the full interval [-Na/2, Nal2\. The unfolded geometry is the one used in our DMRG 
computations, cf. Fig. 1. 
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3. Use of periodic boundary conditions in Fig 2 

In most studies of the 2CK model there is no interaction between the two channels; however, for the results presented in Fig 
2 the two channels are coupled through the periodic boundary conditions. For the remainder of the results presented this is 
not the case. As long as ;r « A/2 this should not affect the results for large enough N since in the limit A —» oo results have 
to be independent of this boundary condition. Obviously, once x approaches A/2 there is no reason to expect that we should 
correctly represent 2CK physics. This can for instance be seen in the inset of Fig. 2(b) where Simp(x,N) does not approach 
ln(2)/2 as x approaches A/2. It is then reasonable to ask why we use periodic boundary conditions since numerically it would be 
simpler to use open boundary conditions. The major reason for this is that the correction to scaling terms appearing in Eq. (9), 
proportional to a, are well understood 1451 . These terms are more correctly seen as arising from ICK physics associated with 
the open boundary in the left panel of Fig. 2(b). In the presence of open boundary conditions several other terms would appear 
dominating the 2CK physics we are trying to study and it is therefore a significant advantage to employ the periodic boundary 
conditions as was done in Fig. 2. Furthermore, we note that Hpsc as dehned in Eq. (2) approach a completely uniform chain in 
the limit 7' —> 1 which is another advantage when extracting S imp. Eor the results presented in Figs. 3,4 the above advantages 
of using HpBc are not pertinent and we have therefore used the more standard Hqbc- 


4. Limiting Behavior of S i„,p as /' —» 1 

As discussed in the main text 5imp(7' ^ 1) = ln(2) - iniguv/gm) using Eq. (4) and (7). This relation takes the view point og 
2CK physics. It is of course also possible to remain with a purely spin-chain view point in which case it follows from Eq. (4) 
and (6) that S imp(7' ~ l,x,N) 2 Ing = ln(2)/2 for 1 <s: x N. These two results just re-express the basic relation 

guvIgiR - (S4) 


discussed in detail in the main text. 


